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ABSTRACT 

A  random  stress-history  \rtiich  is  proportional  to  the  stationary 
response  of  a  single-degree-of- freedom  vibratory  system  to  wide-band 
Gaussian  excitation  is  assigned  a  damage  based  on  the  Palmgren-Miner 
hypothesis  and  an  idealized  S-N  law.  The  damage  acc\amilated  in 
time  T  is  a  random  variable  because  of  the  randomness  in  the  number 
of  "cycles"  and  the  lendomness  in  the  amplitudes  of  the  "cycles." 

The  33iean  and  variance  of  the  damage  are  obtained  by  two  procedures: 
one  \diich  accounts  for  both  sources  of  randomness  and  one  \jhich 
neglects  the  randomness  in  the  number  of  "cycles"  contained  in  the 
interval.  The  two  procedures  give  the  same  aByn5)totic  result  when 
the  bandwidth  shrinks  to  zero.  The  theoretical  results  are  illustrated 
by  curves  computed  for  a  particular  exaniple. 


The  Variance  in  Palmgi«n-Miner  Damage 
Due  to  Random  Vibration 


A  damage  D(T)  can  he  associated  with  an  interval  T  of  a  stationary 
narrow-hand  random  stress-history  s(t)  by  using  the  Palmgren-Miner 
criterion  [^IJ  .  This  damage  is  a  random  variable  taking  on  different 
values  for  each  sample  stress-history.  In  195^  Miles  L  2j  evaluated 
the  expected  value  or  mean  of  the  damage  when  the  stress-history  was  the 
response  of  a  single-degree-of- freedom  vibratory  system  to  ^diite  Gaussian 
excitation  and  the  S-N  diagram  or  fatigue  "law”  for  the  material  was 
assumed  to  have  the  form 


NS'^  =  constant  =  s/' 


We  shall  be  concerned  with  the  same  situation  and  will  evaluate  the  variance 
of  the  damage  in  addition  to  the  mean. 

The  damage  D(T)  when  T  is  large  is  the  siua  of  a  large  number  of 
incremental  damages  each  associated  with  a  single  "cycle'."  The  randomness 
in  D  is  due  to  the  randomness  in  the  amplitudes  of  the  individual  "cycles" 
and  also  due  to  the  randomness  in  the  number  of  "cycles"  contained  in  the 
interval  T.  We  have  made  two  evaluations  of  the  variance;  the  first  takes 
into  accoxint  both  sources  of  randomness  while  the  second  is  an  approximate 
solution  which  considers  only  the  randomness  in  the  amplitudes  and  neglects 
the  randomness  in  the  periods  of  the  "cycles."  The  two  solutions  are  shown 
to  approach  one  another  in  the  limit  as  the  bandwidth  is  decreased  to  zero. 
In  both  cases  the  major  difficulty  in  the  analysis  is  due  to  the  strong 
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correlation  in  the  incremental  damages  of  succeeding  "cycles."  It  has  been 
necessary  to  ruEike  approximations  idiich  are  only  valid  when  T  is  long  in 
comparison  with  the  decay  time  of  the  correlation. 

It  does  not  appear  possible  to  obtain  the  con^ilete  probability 
distribution  of  the  damage  D(T)  for  finite  T  although  the  central  limit 
theorem  can  be  invoked  to  show  that  in  the  limit  as  T-^cothen  the 
distribution  of  D  becomes  normal.  In  this  limiting  sitviation  the  meein  and  the 
variance  are  sufficient  to  completely  characterize  the  distribution  of  D. 

1.  Mean  and  Variance  of  Sums 

Let  the  total  interval  T  be  divided  into  M  equal  subintervals.  With  each 
subinterval  let  an  incremental  damage  be  associated.  For  the  moment 
we  postixjne  the  discussion  of  how  the  incremental  damage  is  to  be  associated. 

The  total  damage 

D  =  Z'  (« 

t-o 

is  the  sum  of  (correlated)  random  variables.  The  mean  and  variance  of  D 
are 

CCD]  2  Eidi] 

L 

.  (3) 

ra^iDl  =  E[Dn-(£[Djf  ^ 

I  J 

Now  since  the  stress-history  is  stationary  the  damage  process  is  also  end  the 
statistical  averages  needed  in  (3)  are  invariant  with  respect  to  a  translation 
of  the  time  axis.  Thus 

E  [d-J  ~  E  LdjJ  ^  E[do2  (4) 

B  [di  dj]  =  Eid^  ^  J 


for  arbitrary  i  sind  j.  The  sums  in  (3)  may  then  be  recast 'as  follows 


^  E  :  D]  =  ,M  ELdJ 

'Mr  f  DJ  =  Ml  ECd,^]  -  cmj/J  (eUJJJ 

Although  the  damages  in  adjacent  subintervals  nay  be  strongly  correlated 
and  become  uncorrelated  when  k  gets  large  enou^  and  thus 

liM  j  rCol^d^']-{E[doJ)^]  -0 

k-^0(O  L  J 

This  will  be  helpful  for  evaluating  (5)  for  large  M. 


2.  Incremental  damage  associated  with  subintervals . 

Ordinarily  the  Palmgren-Mlner  hypothesis  is  used  to  associate  a 
damage 

with  the  i*th  cycle  where  is  the  number  of  cycles  until  fellure  at  the 
constant  stress  amplitude  S^.  as  given  by  the  S-N  diagram  or  by  a  relation 
such  as  (1).  In  order  to  avoid  certain  subtleties  involved  in  determining 
the  peak  amplitudes  of  a  random  process  we  consider  a  slight  modification 

of  the  hypothesis  in  which  we  associate  a  damage 

J 

with  a  zero-crossing  of  the  stress  process  using  the  slope  s  as  a  measure 
of  the  stress-aorplitude .  If  the  expected  frequency  of  the  narrow  band  process 
is  oOo  where 


QCco)  du) 


and  rC'C-')  and  G((-tJ')  are  the  autocorrelation  function  eind  spectral  density 
respectively  of  the  process  (assumed  to  have  zeroi  mean)  then  the  equivalent 
stress  amplitude  we  associate  with  a  zero-crossing  having  slope  a  is 

m 

C^ct 

and  the  niunber  N  in  (8)  is  obtained  from  the  S-K  dlagreun  using  (10)  for 
15.  The  factor  2  in  (8)  arises  because  there  are  twice  as  many  zero-crossings 
as  cycles  in  a  narrow-band  process;  it  could  be  avoided  by  considering  only 
the  zero-crossings  with  positive  elope  but  the  integrations  leading  to  (21) 
are  considerably  eimpler  when  we  associate  damage  with  half-cycles 
rather  than  with  cycles.  The  relation  (3,0)  would  be  strictly  correct  for 
simple  harmonic  motion  at  frequency /;4) .  In  a  narrow-band  process  it 
represents  a  good  approximation  to  the  amplitude  of  the  stress  peak  immediately 
after  (or  before)  the  zero- crossing.  We  believe  that  the  statistics  of  the 
zero-crossing  damage  process  SO’  defined  will  not  differ  slgnificemtly  from 
the  statistics  of  the  peak- associated  damage  process. 

The  two  evaluations  which  follow  are  based  on  two  choices  for  the 
subinterval  duration.  In  the  first  case  the  subinterval  is  taken  to  have 
the  duration  /it  and  eventually  ^t  is  taken  to  approach  zero-  This  procedure 
permits  us  to  take  into  account  the  variation  in  the  periods  of  the  "cycles." 

In  this  case  the  incremental  damage  is  taken  to  be  zero  if  there  is  no  zero- 
crosBlng,  within  the  subinterval  or  to  be  the  value  (8)  if  there  is  a  zero- 
crossing.  In  the  second  case  the  variation  in  the  periods  of  the  "cycles" 
is  neglected  and  the  duration  of  each  subinterval  is  taken  to  be  ly/cJo  ;  i.e._, 
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the  expected  duration  of  a  half-cycle.  The  incremental  damage  is  taltea  to 
be  (8)  but  here  instead  of  using  (10)  we  find  it  more  convenient  to  use 
the  value  of  S  given  by  Rice's  envelope  function  • 

3-  First  ease:  infinitesimal  subintervals. 

The  interval.  T  is  divided  into  equal  subintervals  such  that 


MAi  =  T 


(n) 


We  assume  that  ^it  is  so  small  that  the  stress-history  s(t)  can  be  taken 
as  a  straight  line  throughout  the  interval.  The  fraction  of  samples  which 
wi].l  have  a  zero-crossing  in  a  particular  subinterval  At  can  be  ascertained 
by  considering  the  distribution  of  combinations  8(t)  and  B(t)  where  t  is 
the  time  at  the  beginning  of  the  Interval.  This  distribution  is  described 
by  the  joint  density  function  p(s,s).  Combinations  of  s  and  s  which  involve 
a  zero-crossing  are  those  for  which 

-  5/lt  <  5  <  O  ;  S  >  O 

0  <  S  <  -  SAtj  S  <  O 


(12) 


For  those  samples  having  a  zero-crossing  the  incremental  damage  is  given 
by  (8)j  for  those  without  a  zero-crossing  there  is  no  deiraage.  The  expected 


damage  in  the  subinterval  4t  is  then 

^/V 


ELoiJ  =  )  l^is.s)ds+  ] 

-  60  ^  o  -  SAt 


ist 


Joo 

~  HI  s)  ds 

c>6 


(13) 
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To  evaluate  (13)  we  make  the  assuMption  that  s(t)  is  a  Gauesiaa  process 
so  that  .  ,  ■ 


(11^) 


and  we  assume  that  (1)  is  the  S-N  curve  of  the  material  so  that  using  (10) 


_L 


JAL 

3 i 


(15) 


Inserting  (ih)  and  (15)  into  (13)  yields 


6Jo 


(16) 


Finally  inserting  (l6)  and  (11)  into  (5)  leads  to  the  mean  or  expected 
damage  for  an  interval  T 

E  fD(rj]  =:  EAT  (~)  (17) 

Ox  ' 

where  is  the  expected  number  of  zero-crossings  with  positive  slope;  i.e., 
is  the  expected  or  mean  frequency  la  cycles  per  unit  time 

P/  =  (18) 

eir 

The  result  (17)  although  derived  differently  is  identical  with  that  of  Miles  [2~J  . 
The  above  derivation  can  be  repeated  using  the  square  of  (15)  to  obtain 
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the  mfean  square  damage  associated  vith  the  suhlnterval  At- 


CIS) 


The  derivation  for  the  correlation  terms  E  similar  hut  more  complex. 

The  product  d^dj^  ifl  zero  unless  zero- crossings  occur  in  both  the  zero-th  and 
k-th  Buhintervals  in  which  case  the  product  is  given  by  multiplying  together 
terms  of  the  form  (8)  which  in  turn  depend  (lO)  on  the  slopes  s  at  tl« 
zero- crossings.  The  fraction  of  samples  which  will  have  zero  crossings  in 
both  subintervals  can  be  ascertained  by  considering  the  distribution  of 
combinations  s(t),  s(t),  s{t  +  k  At)  and  8(t  +  kAt)  ^ich  is  described 
by  the  four-dimensional  Joint  probability  density  p(8(t),  6(t),  s(t  +2'), 
a(t  Again  by  Integrating  over  the  subregion  where  both  subintervals 

have  zero-crossings  we  find,  analogous  to  (13) 

o<T 

F  Zdo 4]  =  Mi /  j  \  i4  (13) 

C?C>  —  <?<} 


for  k'/O.  To  evalxiate  (19)  we  nsake  the  a&sumption  that  B(t)  is  a  Gaussian 
process  so  that  j  3J 


exj} [-  ^  (Aa^ 


..a 


'df  yj  ( 90) 


where  the  yl- parameters  depend  on  the  autocorrelation  of  the  stresB  process  and 
are  evaluated  in  Sec.  4-.  Again  we  assume  that  the  jEncxenental  damages  are 
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given  "by  (15)  but  here  to  render  the  integral  tractable  we  nake  the  additional 
assumption  that  the  exponent  cx  is  an  odd  positive  Integer.  In  this  way  we 
find  [Uj 


n  2. 


A 

A 


(X+  2  /  ( 

22. 


PA, 


■  L  ■ 
2^ 


A 


id- 


A 


(21) 


iL 


Finally  to  obtain  the  variance  in  the  total  damage  we  Insert  ( 21) ^  ( 18) 
and  (16)  into  (5).  At  this  time  we  also  let^t-r^-O  thereby  converting  the 
siujanation  into  a  Riemann  integral.  With  kylt  =  'Z'  we  find 


Further  evaluation  recLuires  a  specific  choice  for  the  autocorrelation  fhnetioa 
R(  2^)  in  order  to  specify  the  -  parameters  and  the  hypergeometric  function. 

k.  Specialization  to  the  response  of  a  single-degree~of-feedom  system. 

We  limit  our  discussion  to  the  case  \rtiere  the  stress  history  8(t)  le 
proportional  to  the  response  of  a  lightly-damped  single-deg3?ee-of- freedom 
oscillator  when  excited  by  stationary  white  noise;  i.e.^  s{t)  is  taken  to 
satisfy  the  differential  equation 

5  -t  2.2^ cOy,  5  S  =  -fc-t)  (23) 


where  is  the  damping  ratio  and  u},,  is  the  undamped  natural  frequency.  When 
the  excitation  f(t)  is  stationary  idiite  noise  the  autocorrelation  function  of 


dz 
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the  response  is 

R(r;  ^  sm/oz)  ^  (2^) 


where  p  =  /l  - ^  6J„  is  the  damped  natural  frequency  and  CI^  ^  is  the  mean 
square  stress.  The  parameters  are  obtained  from  R(7^)  according  to  the 
following  definitions 


A  =•  L  R'o]  -  RlV^J\_  R"(o)^~  R"(T)^J-t£Rcz)  lRco)R’(o)-R(T)R"('Zyi}  R(T) 


Aaa'-'  -  -  RA)  C  R(o/-  /RCTJ~]  -  R(o)R'[zf 

R"'tr)iRic)A  Acz)'2  Rn)R'cz)^ 


Using  (24)  we  find  that  the  expected  frequency  of  the  narrow  hand  process^ 
of  (9)^  is  Just  the  natural  frequency  o)^  .  It  is  thus  possible  in 
principle  to  insert  (25)  and  (24)  in  (22)  and  evaluate  the  variance  of  the 
damage.  The  integration  appears  however  to  involve  formidable  difficulties. 
We  therefore  discontinue  the  exact  evaluation  and  consider  an  approximation 
for  small  damping.  We  use  the  following  small  ^  approximations  to  (25)  for 
use  in  (22) 


/i 
i.  1 


) 


Alo:  (/- 


(T 


(  1-  e 


2.  ^  J  1  i- 


(26) 


A,/ 


S' 


ucJlj  ~  ccl  L 


10- 


The  argument  of  the  hypergeometric  function  thus  oscillates  with  period 
u),-  .  The  integrand  in  (22)  has  peaks  at  and  fal].8  to 

zero  at  (2n  +  1)'!t/2u)^,  ,  n  =  1,2^...  .  iXirthenmore  the  integrand  is 

negligible  in  the  range  0  <7r<  ?r/26<)n.  We  have  been  able  to  obtain  a 
good  approximation  to  the  Integral  by  substituting  a  smoothed  integrand 
in  which  the  expression  in  braces  in  (22)  is  replaced  by 

[  F(-f.-fs  -  ij  (S7) 

and  the  limits  of  integration  are  changed  to  vcj2Q^  and  i-O  .  The 
approximations  involved  are  good  for  1  and  ^  »  1.  The  result 

is 


(28) 


where  the  f- quantities  are  given  by 


CO 


'■  cj„dx 

‘X  i  ^  j 

c  j  ^  A  Z 

(29) 

4  ~ 

1  (l  -/-  X  ) 

1  tr.  1  0-1-  X/2  ) 

and  are  tabulated  in  Table  I .  A  partial  check  on  the  accuracy  of  approximation 
was  made  by  comparing  the  smoothed  integral  used  above  with  the  results  of  a 
numerical  integration  of  the  exact  integral.  For  (X  =  9  and  /'■  =  l/60  the 
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discrepancy  vas  less  than  one  part  in  300.  The  errors  tend  to  increase 
rapidly  wither  “but  for  the  range  tabulated  in  Table  I  the  result  (g6)  is 
probably  satisfactory  for  engineeriixg  purposes  for  O.05  (Q>10)  if  the 
expected  number  of  cycles  a^^t  is  very  large  compared  with  Q  =  l/2^  >  In 
this  range  very  little  additional  error  is  made  by  neglecting  the  fg  and  f^ 
terms  in  comparison  with  the  f^^  term  in  (S8) .  Thus  in  the  range  indicated 
we  may  use 


which  indicates  a  linear  growth  of  the  variance  with  T  and  an  inverse 
dependence  on  the  damping  ratio ^  for  a  given  material  and  a  fixed  mean 
square  stress  level. 

If  we  denote  the  standard  deviation  of  the  total  damage  by  we 
can  combine  the  results  (17)  and  (30)  for  the  wssn.  and  variance  into  the 
following  ratio 


which  indicates  that  the  relative  variance  for  a  given  material  decreases 
in  inverse  proportion  to  the  square  root  of  the  product  of  T  and  the  bandwidth 

2  ylf  Ol^Q). 

5.  Second  case;  half-cycle  subintervals • 

Here  we  consider  the  stress  process  B(t)  to  consist  of  a  sequence  of 
half-cycles  each  of  duration  'iT/cO^  ;i.e.,  we  neglect  the  random  variation  in 
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periods.  With  each  half-cycle  the  incremental  damage  (8)  is  l/2N  vhere 
the  S-N  diagram  or  (1)  is  used  to  relate  N  to  a  stress  amplitude  S  associated 
with  the  half-cycle.  For  this  purpose  we  find  it  convenient  to  use  Rice's 
envelope  function  S(t) .  ^or  a  narrow-hand  process  Rice  (  3  f  has  shown  that  the 
first  and  second  order  prohahility  densities  for  S(t)  are 


i  h  2 

where  is  the  mean  square  of  the  process  s(t)  and  the  quantities  A  =  -B 

and  B  are  functions  of  the  time  interval  t^  -  t^  and  depend  on  the  spectral 

density  of  s(t)  .  They  will  he  described  later. 


The  expected  value  of  the  incremental  damage  is 


'c.  i  4, 


^  LS)  dS 


(33  j 


Substituting  from  (1)  and  (32)  yields 


I 

V  S', 


r  oc/^) 


Note  that  this  is  equivalent  tp  (16)  with  At  =  TT/cJo  •  This  is  essentially 

the  teclmique  used  by  Miles  ^2  |  in  deriving  (IT)  •  An  exactly  similar  evaluation 
2 

using  (l/a?)  in  place  of  l/2N  in  (33’)  leads  to 


E  ['.i'-'J  =  Z  ( 


I  ^  C  /  ) 


(35) 
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which  is  equivalent  to  (l8)  with  t  ^iTfcOo  • 

Turning  next  to  the  correlation  terms  in  (5)  we  have 

Ez  11‘^J'k  ]  -  \  ^  ^  IP^  S.;  Sj^)  ^ Sa  O/^  (36)' 

which  on  substitution  from  (l)  and  (32)  yields 


E'  L  (Po  = 


^  -  1'  H 


(37) 


which  is  now  q.ulte  different  from  (21)  because  (37)  is  the  correlation 
between  two  half-cycles  separated  by  kT/cJ^  while  (21)  is  the  correlation 
between  two  infinitesimal  intervals  separated  by  k/).t.  Finally  to  obtain 
the  variance  we  insert  (37)),  (35))  and  (3^)  into  (5)  to  obtain  the  variance 
after  M  half-cycles  (T  =  MiT'/^o  =  M  i^J/2) 


Var  [  D(T;J 


P  O' 


Further  evaluation  requires  a  specific  choice  for''s^ctral  density  G(<^J) 
of  the  process  in  order  to  specify  the  parameter  B  and  hence  the  hypergeometric 


function. 


6 .  Special! aatlon  to  the  response  of  a  slngle-degree-of-feedom  systeia  . 

We  limit  our  dlBcitisBiion  again  to  the  process  deserihed  in  Sec.  h.  The 
spectral  density  of  the  process  is 


(co  '  -  '  O,.!  ) 


-r 


(39) 


and  the  parameter  B  follows  from  the  following  operation 

^2-  1-  ,2. 

13  -  y^\3 


^  c<? 

)  QiO-o)co$  ctK)-u)r,)T  (^<-0 

‘'-4 

Q  (Caj)  51*]  (I'O-  6(1^ I  X  d<A) 


c 


(4o ) 


This  is  quite  complicated  in  general^  hut  here  we  are  only  interested  in 
times  '>■  of  the  form  kir/uJ,.,  and  moreover  we  will  again  accept  the  same  type 
of  small  damping  approximations  used  In  Sec.  4.  Under  these  circumstances 
we  find  simply 


(?+l) 


and  the  term  in  the  traces  heing  summed  in  (38)  becomes 


!  fX 

r  '5: 


.  -  ;r  ^  i  j 


(42) 


At  this  stage  there  is  considerable  similarity  between  the  term  being  summed  (4.2) 
and.  the  smoothed  integrand  term  (2?)  .  In  fact  the  sum  may  be  considered  as  a 
crude  attempt  to  approximate  the  integral.  The  hypergeometric  functions  in  (  38) 
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can  be  expanded  In  series  and  the  order  of  summtions  interchanged.  The 
k- summations  involve  simple  geometric  progressions  which  can  be  summed.  When 
M  is  large  and  .^is  small  the  resulting  series  can  be  recognized  as  equivalent 
to  the  functions  defined  in  (29)“  In  this  way  we  obtain 


V - 


-  -h  f 

J 


(*+3) 


for  the  variance  when  ^  is  small  and  the  expected  number  of  cycles  T  is 
large.  Again  if  is  large  (even  though  cr^  is  small)  we  can  dispense 

with  the  final  two  terms  and  use  the  asymptotic  form 


yar[Dirj] 


which  should  be  compared  with  (30 ) .  We  note  that  the  discrepancy  between  the 
first  case  which  accounted  for  the  random  variations  in  periods  and  the  second 
case  which  neglected  this  source  of  variance  is  small  for  light  damping  and  that 
the  ratio  of  the  two  expressions  (30)  and  (M:-)  approaches  \mity  as  ^->>0. 

7 .  Example 

Consider  the  system  shown  in  Fig.  1  in  which  the  vehicle  has  a  stationary 
random  acceleration  with  uniform  spectral  density  0.5  g^/cps.  With  the 
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following  data 

4.0" 

0.25" 

10.3  X  lO^psi 

6.09 

5 

2  X  10  psl 

7.28  X  10^ 

lb  ^ec^/in 

tbe  natural  frequency  of  the  system  is 

=  465  rad/sec  (73-9  cps) 


length  of  cantilever, 
h,  side  of  square  cross  section 
E,  75S-t6  aluminum  alloy 

^  II  II  II 

\ 

tl  11  II 

m  p  mass  of  one  cubic  inch  of  steel 


and  the  rms  stress  level  of  the  narrow-beind  stress  response  in  the  extreme 
fibers  at  the  root  of  the  cantilever  is 


Z3'dO 


psi 


(^5) 


where ^ is  the  damping  ratio  of  the  system.  With  these  values  it  is  now  poesible 
tOy^the  expected  damage  (17)  and  the  standard  deviation  of  the  damage  (31) 
as  functions  of  ^  and  T-  The  value  2-32  for  fq(P^)  i^  obtained  from 
interpolation  in  Table  I.  The  results  are  shown  in  Fig.  2  for  four  different 
values  of  system  damping.  The  expected  damages  appear  as  straight  lines  with 
unit  slopes  Indicating  linear  growth  with  time.  The  variance  in  the  accumulated 
damage  is  suggested  by  the  curves  showing  E  [jdJ  +  (j^  .  Note  here  that  the 
primary  effect  of  the  dandling  is  through  its  action  in  setting  the  rms  stress 
level  ((455)  while  a  secondary  effect  is  its  action  in  controlling  the  deviation  of 
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the  damage  through  the  correlation  of  the  damages  of  successive  cycles  (  3I)  • 

It  is  also  of  interest  to  compare  "the  two  solutions  (30)  and  ('44)  for  tht, 
variance  in  the  range  covered  by  Fig.  2.  We  find  a  2%  discrepstncy  in  the 
variance  for  Q  =  10  and  2- 5^  for Q  ==  100.  The  discrepancies  in  the  deviation 
('7,3  would  be  about  half  of  these  values. 

The  Palmgren-Mlner  criterion  for  failure  is  D  ss  1.  Fig.  2  indicates 
that  when  the  average  damage  reaches  unity  there  is  actually  a  distribution 
of  damage  across  the  ensemble  of  sample  histories.  There  will  correspondingly 
be  a  distribution  of  time- to -failure  Tp.  This  distribution  is  unknown  but 
the  central  limit  theorem  can  be  invoked  to  show  that  It  also  becomes 
asymptotically  Gaussian  as  T  ->'K)  and  that  the  plus  and  minus  one-signui 
limits  for  Tp  are  asyn^ptotically  the  points  where  the  plus  and  minus  one- sigma 
curves  for  damage  cross  D  =  1. 
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Tatle  1  Quantities  in  Equation  (28)  Which  Depend  on 
the  Fatigue  Law  Exponent  r  \  . 


1 

0 .04l4 

0.00323 

0 .0796 

3 

0.369 

0.0290 

0.212 

5 

1.280 

0.0904 

0.679 

7 

3.72 

0.223 

2.33 

9 

10.7 

0.518 

8.28 

11 

31.5 

1.230 

30.0 

13 

96.7 

3.06 

111.2 

15 

308. 

8.11 

415. 
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CAPTIOMS  FOR  FIGURES 


Fig.  1  (a)  Random  stress  history  at  A  is  due  to  random  vibration  of 

vehicle,  (b)  Schematic  excitation  -  response  diagram. 

Fig.  2  Pa Imgren- Miner  dauiage  at  root  of  cantilever  beam.  The  mean 

damage  expected  is  shown  together  with  the  plus  and  minus 
oxie-sigm  limits. 
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s(t)  = 
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